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Abstract 



We consider (for the first time) the ratios of doubly heavy baryon masses (spin 3/2 over spin 1/2 and SU(3) mass-splittings) using 
double ratios of sum rules (DRSR), which are more accurate than the usual simple ratios often used in the literature for getting 
the hadron masses. In general, our results agree and compete in precision with potential model predictions. In our approach, the 
a.s corrections induced by the anomalous dimensions of the correlators are the main sources of the Hgg - Sqq mass-splittings, 
I which seem to indicate a 1/Mq behaviour and can only allow the electromagnetic decay Sgg — > Sqq + y but not to 'B.qq + n. Our 
results also show that the SU(3) mass-splittings are (almost) independent of the spin of the baryons and behave approximately like 
'1 IMq, which could be understood from the QCD expressions of the corresponding two-point correlator. Our results can improved 
O including radiative corrections to the SU(3) breaking terms and can be tested, in the near future, at Tevatron and LHCb. 
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1. Introduction 

(—1 

Q_|In a previous paper IH], we have considered, using double ra- 
' ItiosilSl of QCD spectral sum rules (QSSR) QS] (DRSR), 
^^the splittings due to SU(3) breakings of the baryons made with 
f~| pne heavy quark. In this paper, we pursue this project in the 
'""'case of doubly heavy baryons. The absolute values of the dou- 
bly heavy baryon masses of spin 1/2 (Hqq = QQu) and spin 
^ 3/2 (Sgg = QQu) have been obtained using QCD spectral sum 
rules (QSSR) (for the first time) in lH] with the results in GeV: 
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3.58(5) . 
3.48(6) , 



•and in IH: 



6.86(28) . 



: 10.33(1.09), 
9.94(91), (1) 



(2) 



More recently ifTH [12I1 . some results have been obtained us- 
ing some particular choices of the interpolating currents. The 
predictions for M-^-^, and M=^_. are in good agreement with the 
experimental candidate Me,, = 3518.9 |13]. We shall also im- 
prove these previous predictions by working with the DRSR 
for estimating the mass ratio of the 3/2 over the 1/2 baryons 
and shall compare them with some potential model predictions 
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2. The interpolating currents and the two-point correlator 

For the spin 1/2 QQq baryons , and following Ref. |9], we 
work with the lowest dimension currents: 



J^a = [(QlCysqp) + biQlCqp)y5] Qa- 



(3) 



where q = d,s are light quark fields, Q = c,b are heavy quark 
fields, b IS, a priori an arbitrary mixing parameter. Using the 
/7-stability criterion of the QSSR results for the masses and cou- 
plings, the optimal values of these observables have been found 
for: 

^ = -1/5, (4) 

in the case of light baryons UtIi and in the range ifH [l"8l - l20ll : 

-0.5<fe<0.5, (5) 

for non-strange heavy baryons . The corresponding two-point 
correlator reads: 



S{q) 



i d^xe^^'-'' {0\TJ.Jx)J.JOm 



= qFi+F2, (6) 
whereH: q = t^. The invariants Fj ( j - 1 , 2) obey the dispersion 



relation: 



Fj(q^) 



f 

J(21\ 



dt 



-lmFj(t) + . 



hlMQ+mqf t - q - le 71 

where . . . indicate subtraction constants and -q^ = Q^ > 0. 
For the spin 3/2 QQq baryons, we also follow Ref. ^ 
work with the interpolating currents: 



(7) 



and 



4.^ = ^|\e,,p,[2{Qlcy'dp)QA + {QlcrQp)qA\ (8) 



We use the notation in the Landau and Lifchitz's book. 
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The corresponding two-point correlator reads: 

J ~0 -Q 



^'"'(§^1+^2) + ... , 



3. The two-point correlator in QCD 



(9) 



The expressions of the two-point correlator using the previ- 
ous interpolating currents have been obtained in the chiral limit 
ma - and including the mixed condensate contributions by 
In this paper, we extend these results by including the lin- 
ear strange quark mass corrections to the perturbative and (ss) 
condensate contributions. We shall use the same normalizations 
as in d. 

For the spin 1/2 baryons, these corrections read: 
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For the spin 3/2 baryons, these corrections read: 
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with: 



(11) 



(12) 



4. Form of the sum rules and QCD inputs 

We shall work with the exponential sum rules ifTi I24II25I1 : 



Tiir) = I dt e-" -\mFi{i) , (/ = 1, 2) 
n 



(13) 



(r =\ jM^ is the sum rule variable) from which, one can derive 
the following ratios: 



/■CO 

( dt t g-'^ ImF,(0 
r dt e-"\mFi{t) 
r dt e-" ImFtit) 



(/= 1,2), 



(14) 



used in the sum rule literature for extracting the baryon masses. 
We parametrize the spectral function using the standard duality 
ansatz: "one resonance" H- "QCD continuum". The QCD con- 
tinuum starts from a threshold tc and comes from the disconti- 
nuity of the QCD diagrams, which is consistent with a matching 
of the QCD and the experimental sides of the sum rules for large 
t. The value of t^ is not arbitrary as its value obtained inside the 
region of fc-stablity of the sum rule is correlated to the ground 
state mass and coupling f2?]. This simple duality model has 
been successfully tested in the literature when a complete data 
for the spectral functions are available, like e.g., e^e' to 1=1 
hadrons or to charmonium data [8]0. Transferring the QCD 
continuum contribution to the QCD side of the sum rules, one 
obtains the finite energy inverse Laplace sum rules: 



dt e-" -ImF2(0 . 
n 



dte'" -lmFi(f) , 

71 



(15) 



(10) where /t^w and M^w are the heavy baryon residue and mass 
from which one can derive the FESR analogue of the ratios of 
sum rules. Consistently, we also take into account the SU(3) 
breaking at the continuum thresholdQ: 



sua) 



- ( Vf^lsf/(2) = yfC) 



(16) 



fhs is the running strange quark mass. As we do an expansion 
in nis, we take the threshold tg = for consistency. niQ is the 
heavy quark mass, which we shall take in the range covered by 



^More involved parametrizadons of the continuum can also be proposed 
(see e.g. (22l for non-resonant final states within ChPT or (23ll for a t-dependent 
tc model.). However, it is easy to check in the harmonic oscillator model dis- 
cussed in 1 23] that the 5% uncertainties induced e.g. by a ^-dependent contin- 
uum model on the ratio of moments will be negligible in the double ratio of 
sum rule (DRSR) defined in Eq. i20\ as the leading corrections coming from 
light-quark flavour-independent terms will largely cancel out. This cancella- 
tion of the QCD continuum contribution will be signaled by the large range of 
(c-stability region obtained in our analysis. 

'^As we have done an expansion in terms of mj, the quark threshold has to 
be taken at 4m^ but not at (2mg + f/ij)^. 



the running and on-shell mass (see Table [T]) because of its am- 
biguous definition when working to LO. At the r-stabihty point 
of the FESR analogue of the ratios of sum rules, one obtains: 



0'= 1,2) 



(17) 



These predictions lead to a typical uncertainty of 10-15% jol, 
which are not competitive compared with predictions 
from some other approaches, especially from potential models 
131 ■ In order to improve the QSSR predictions, we work with 
the double ratios of finite energy sum rules (DRSR)B 




1,2) 



^sd ^ '221 



(18) 



which take directly into account the SU(3) breaking effects. 
These quantities are obviously less sensitive to the choice of the 
heavy quark masses and to the value of the continuum threshold 
than the simple ratios % and iRaiEl For the numerical analysis 
we shall introduce the RGI quantities fi and rhg ll26ll : 



{qq}(T) = 
{qGq){T) = 



(-logV^Ap 

-/i^,(-logV?Ap 

-A^X-logV7A)''"^' 



Mi 



(19) 



where fix - -(1/2)(11 - 2n/3) is the first coefficient of the 
j6 function for n flavours. We have used the quark mass and 
condensate anomalous dimensions reviewed in |8]. We shall 
use the QCD parameters in Table [1] 

- We shall not include the Ijq^ term discussed in |j27 , 28 ] ,which 
is consistent with the LO approximation used here as the latter 
has been motivated for a phenomenological parametrization of 
the larger order terms of the QCD series. 

- We have used the value of k = {ss)/{dd} from |1] which 
we consider as improvements of the ones from light meson 
systems id, 2^31|, where the one from the scalar channel 
suffers from the unknown nature of the k meson, while the 
one from the pseudoscalar channel depends on the theoretical 
appreciation of the n' meson contribution into the spectral 
function l^ll^]. However, the different estimates agree each 
others within the errors. To be conservative, we have multiplied 
the original error in |[ll] by 2. 

- For the gluon condensate, we have used the estimate from 



heavy quarkonia and e^e" data 112 iL l24i l32l - l37ll . We do 
not expect that the estimate from r-decays is reliable as its 
contribution acquires an extra-o-s term in the t width compared 



^Analogous DRSR quantities have been used successfully (for the first time) 
in for studying the mass ratio of the 0+^/0"* and 1^+/1 B-mesons, in @l 
for extracting fsj fs, in 15] for estimating the D ^ K/D ^ n semi-leptonic 
form factors and in |6] for extracting the strange quark mass from the e^e" — » 
/= 1,0 data. 

^One may also work with the double ratio of moments M„ based on differ- 
ent derivatives at = |3]. However, in this case the OPE is expressed as an 
expansion in I/ihq, which for a LO expression of the QCD coiTclator is more 
affected by the definition of the heavy quark mass to be used. 



to the one in the two-point correlator fss'], while its value, 
in this process, can also be affected by the treatments of the 
large order PT series |27, [32I1 . However, the effect of the 
gluon condensate is not important in our analysis of the SU(3) 
breaking as it disappears like some other flavour-independent 
contributions in the DRSR. 

- For the heavy quark masses, we use the range spanned by the 
running MS mass mQ(MQ) and the on-shell mass from QSSR 
compiled in page 602, 603 of the book in fl- 



Table 1 : QCD in put p arameters. The values of A, m j and /j^ have been obtained from 
a,(Mr) = 0.325(8) [3J and from the mnning masses: m,(2) = 96.1(4.8) MeV and m,,(2) = 
5.1(2) MeV|23]. The original enors for k and {a^G^} have been multiplied by 2. 



Parameters 



Values 



Ref. 



A{nf = 4) 


(324+ 15)MeV 


[13, 32] 


A{nf = 5) 


(194± 10)MeV 


U3,32] 


mj 


(114.5 ± 20.8) MeV 


[6,l,i3,2g] 


fid 


(263 ± 7) MeV 


[8,22] 


K = (ss}/(dd} 


(0.74 ± 0.06) 


Ll,8,29,3r] 


Ml 


(0.8±0.1)GeV2 


[3, r7,32] 




(6 ± 2) X 10-2 GeV* 


[21, 24, 32-37] 


nic 


(1.26 ~ 1.47) GeV 


[a,J3,29,37,40] 


nib 


(4.22 ~ 4.72) GeV 


Il,i3,29,37,40] 



5. The E^g/Sgg mass ratio 
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Figure 1: Charm quark: fo-behaviour of the dilferent DRSR given r = 0.8 



GeV"2 and t,. 



(green); r^^' dashed line (blue). We have used ntc = 1.26 GeV and the other 
QCD parameters in Tabled 



25 GeV^. 



dot-dashed line (red); 



continuous line 



We extract the mass ratio using the DRSR analogue of the one 
in Eq. ( fTSl ) which we denote by: 



3/1 



1,2 



3/1 _ '22}_ 
''21 - ^1 
'^21 



(20) 



where the upper indices 3 and 1 correspond respectively to the 
spin 3/2 and 1/2 channels. We use the QCD expressions of the 
two-point correlators given by |9] which we have checked. We 
notice like |9] that the mixed quark condensate contribution has 
a term which behaves like 1 / (where v is the heavy quark ve- 
locity), which signals a coulombic correction and would require 
a complete treatment of the non-relativistic coulombic correc- 
tions which is beyond the aim of this paper Therefore, in our 
analysis, we truncate the QCD series at the dimension-4 con- 
densates until which we have calculated the m, corrections. We 



3 



shall only include the effect of the mixed condensate (if nec- 
essary) for controlling the accuracy of the approach or for im- 
proving the T or/and tc stability of the analysis. 

The charm quark channel to lowest order in 

Fixing T = 0.8 GeV"^ and tc - 25 GeV^, which are inside 
the T- and fc-stability regions (see Figs. |2]and|3]l, we show in 



3/1 



Fig. [T]the Z7-behaviour of r"*^' which shows that r^ and Tj' ' are 



very stable but not r^^^ ■ We then disfavour r^^- Some common 
solutions are obtained for: 



.3/1 



b ^ -0.35 , and +Q2 , 



(21) 



which are inside the range given in Eq. (|5]l. For definiteness, 
we fix = -0.35 (the other value b = 0.2 gives the same result) 
and study the r-dependence of the result in Fig. |2] We have 
checked in Fig. |2]3 that the inclusion of the mixed condensate 

3/1 

contribution does not affect the result from r- (i — 1,2) ob- 
tained by retaining only the dimension-4 condensates (Fig. |2^) 

^/ 1 

but affects the one from r^l^ . Therefore, we shall only retain the 
results from r ^' (/ =1,2) and show their fc-dependence in Fig. 
[3] The large stability in tc confirms our expectation of the weak 
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Figure 2: Charm quark: a) r-behaviour of Tj'': dot-dashed line (red), rj'' : 
continuous line (green) and Tp' dashed line (blue) with h = -0.35 and = 25 
GeV^. b)the same as a) but when the mixed condensate is included. 
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Figure 3: Charm quark: /r-behaviour of rj'': dot-dashed line (red) and r,'' : 
continuous Une (green) with b = -0.35 and r = 0.8 GeV"-. 

fc-dependence of the DRSR and then on the non-sensitivity of 
the results on the exact form of the QCD continuum including 
an eventual slight f-dependence of tc advocated in |23]. In these 
figures, we have used nic - 1.26 GeV. We have also checked 
that the results are insensitve to the change of the charm mass 



to me = 1 .47 GeV. From these previous analysis, we deduce to 

3/1 

lowest order from r. ' (; = 1,2): 



M: 



- = 0.9994(3) . 



(22) 



The tiny error is the quadratic sum due to {osG^}, nic and a^. 

The bottom quark channel to lowest order in as 
We extend the analysis to the case of the bottom quark. 




Figure 4: Bottom quark: a) r-behaviour of rj'': dot-dashed line (red), r^'': 
continuous line (green)and r^j' : dashed line (blue) with b = -0.35 and fj. = 100 
GeV~; b) the same as a) but when the mixed condensate is included into the 
OPE. 



1.00001 f - 
1.00001 

1-, 

1 \ 

0.999997 ' 

\ 

0.999995 , 

0-999992 V . ^ ■ 

80 100 120 140 160 180 200 

Figure 5: Bottom quark: -behaviour of rj'': dot-dashed line (red) and r^'' : 
continuous line (green) with b = -0.35 and r = 0.6 GeV"^. 

The corresponding curves are qualitatively similar to the charm 
quark one. We take b - -0.35 like in the case of the charm 
quark. The r-stability is reached for t > 0.6 GeV"^ as shown in 
Fig. 3] where we also see that r^^^ is more affected by the mixed 

3/1 

condensate contributions than r. ' . Therefore, we shall elimi- 
nate it from our choice. Another argument raised later about 
the radiative corrections does not also favour rjj'. In Fig. |4l 

3/1 

we study the fc-stability of r. which is reached for tc > 95 

9 3/1 

GeV . Within these optimal conditions, one deduces from r- 
to lowest order: 



'-bh 

Me.. 



1.0000. 



(23) 
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Estimate of the Oia^) corrections 

Radiative corrections due to Os are known to be large in the 
baryon two-point correlators BlTlklll . However, one can easily 
inspect that in the simple ratios "R^ and !R' these huge correc- 
tions cancel out, while its only remain the one induced by the 
anomalous dimension of the baryon operators. Including the 
anomalous dimension 7= 2(resp -2/3) for the spin 1/2 (resp 3/2) 
baryons [17], one can generically write the PT expressions of 
the moment sum rule defined in Eq. ( fT3l l to leading order in 
t/rn^, which is a crude approximation but very informative: 



pert 



(asir)) A Ai T- 



(24) 



where /3i is the first coefficient of the /3-function; A, is a known 
LO expression; Kj is the radiative correction which is known 
in some cases of light and heavy baryons [[l7, 41]. From the 
previous expression in Eq. (l24l i. one can derive the ratio of sum 
rules defined in Eq. (O and then the DRSR in Eq. (fTsT i: 
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Figure 6: Charm quark: a) r-behaviour of r'^'^: dot-da.shed line (red) and 
r^'' : continuous line (green) with h = -0.35, = 25 GeV^ where radiative 
corrections have been included. Bottom quark: b) the same as in a) but for the 
bottom quark. We use h = -0.35 and = 100 GeV- 

3/1 

It is important to notice for r. ' that the radiative correction has 
been only induced by the ones due to the anomalous dimen- 
sions, while the one due to Kj cancels out to this order This is 

3/1 

not the case of where the radiative correction is only due to 
K2 - Ki and needs to be evaluated which is beyond the aim of 
this letter Therefore, in the following, we shall only consider 
the results from r^^'. In our numerical analysis, we shall in- 
clude the as correction into the complete LO expressions of the 
correlators. We show the r-dependence of the DRSR in Fig. |6] 
We shall take the range of r-values where the LO expressions 
have T-stability which is (0.7-1) GeV"^ for charm and (0.5-0.8) 
GeV"^ for bottom (see Figs. |2]and|4]i. One can also notice that 
the NLO DRSR for charm presents a r-extremum in the above 
range (0.7-1) GeV"^ of t rendering its prediction more reliable 
than for the bottom channel case. We can deduce : 

1.0167(10)„/16)„,,., ^ = 1.0019(3)„,(2)m,.(26) 



This would correspond to the mass-splittings (in units of MeV): 

Mh-,. - Mh„. = 59(7) , Mhj,^ - Ms,, = 19(3) , (27) 

if one uses the experimental value 3.52 GeV of the Ecc mass 
which agrees with the QSSR prediction in Eq. ([T]). For the Eht, 
mass, we have used the central value 9.94 GeV in Eq. ^ . The 
ccq mass-splitting is comparable with the one of about 70 MeV 
from potential models 1 10., ,14.1 but larger than the one of about 
24 MeV obtained in 01611 . The bbq mass-splitting also agrees 
with potential models and seems to indicate a 1 /Mh behaviour 
which is also seen on the lattice |42]. Our result excludes the 
possibility that M-^'^^ > M-^qq + m„, indicating that it can only 
decay electromagnetically: 



27- 



(28) 



A future discovery of the S*^ and HJJ^ can infirm or support our 
predictions given to that order of QCD perturbative series. We 
consider the previous results as an improvement of the former 
ones deduced from the mass values in Eq. ([TJ obtained by |@]: 



1.03 + 0.03 



Ms,,, 



1.04 ±0.23 . 



(29) 



6. The figg/Sgg mass ratio 
We use the DRSR in Eq. ( fTSl l where their QCD expressions 
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Figure 7: ilcc/^cc '• a) r-behaviour of r|''(cc): continuous Une (green) and 
rj'^(cc): dot-dashed Hne (red) in the charm quark channel for h = -0.35, = 12 
GeV- and = 1.26 GeV. b) /^-behaviour of rj''(cc) for t = 1 GeV"^: dot- 
dashed line (red) 

can be obtained from the one of the two-point correlator in 
and the new quark mass corrections in Eq. ( fTOb . One can also 
deduce from Eq. jTM that the light-flavour independent radia- 
tive corrections including the one due to the anomalous dimen- 
sions disappear in the SU(3) breaking DRSR, while the most 
relevant radiative corrections are the one corresponding to the 
nis and {ss} terms which are beyond the scope of the LO anal- 
ysis in this paper We show in Fig. |7^ the r-behaviour of the 
DRSR for Mc = 1.26 GeV and b = -0.35 for a given t^ = 10 
GeV^. We have not shown rj2(cc) which is the lesser stable 
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Figure 8: iljj /S^j : a) r-behaviour of r^'{hb): dot-dashed line (red) in tlie 
-0.35, = 100 GeV^ and mj, = 4.22 GeV. b) 
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-behaviour of r'fQyb) for t = 0.5 GeV ^: dot-dashed line (green) 



Figure 9: £1*^/2*^ : a) r-behaviour of rj''(cc)*: dot-dashed line (red) and 
r^^(cc)': continuous line (green) in the charm quark channel for tc = 20 GeV'^ 
and nic = 1.26 GeV. b) f^-behaviour of rj''(cc)* and r2''(cc)*for r = 0.7 GeV"^ 



among the three. We see that the most stable resuh is given by 
r^y{cc). We show in Fig. [TJ) the ff-behaviour of r'^^icc) for a 
given T = 1 GeV"^. We deduce from the previous analysis: 



rf(cc) ^ 



= 1.026(5)„,,.(2).,,(4)„ 



(30) 



where the sub-indices indicate the different sources of errors 
(the parameters not mentioned induce negligible errors). This 
ratio corresponds to : 



Mn„ - Me„. = 92(24) MeV 



(31) 



where we have taken the experimental value M-^^ ^ 3.52 GeV 
from fl^l ■ The errors induced by the other parameters in Table 
[T] are negligible. We perform an analogous analysis in the b- 
channel, which we show in Fig. [8] In this case, we obtain: 

r'\bb) - 1.0049(7)^,(3),-,(10)„, , (32) 

which corresponds to: 

Ma„„ - Me,,,, = 49(13) MeV , (33) 

when we take the value M^,,,, - 9.94 GeV from Our results 
indicate an approximate decrease like 1/wq of the mass split- 
tings from the c to the b quark channels. This behaviour can 
be qualitatively understood from the QCD expressions of the 
corresponding correlator, where the OTj corrections enter like 
Wj/wg, and which can be checked using some alternative meth- 
ods. 



7. The iigg/3gg mass ratio 

We pursue our analysis for the spin 3/2 baryons. The QCD ex- 
pression of the ratios of moments can be obtained from the ones 
of the two-point correlator in |9] and the new mass corrections 
given in Eq. (fTTT i. Including the contributions of the dimension- 
4 condensates, we show your analysis in Fig. |9] One can see in 
Fig. |9^ that r?'' and are quite stable versus t from t > 0.4 




0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



100 
sd. 



Figure 10: 42* la* : a) r-behaviour of r". (bb)': dot-dashed line (red) and of 

bb bb A 

rf^(bc): continuous hne (green) in the bottom quark chaimel for t,- = 100 GeV^ 
and mb = 4.22 GeV. b) Jr-behaviour of rfibby for r = 0.5 GeV"- 



GeV ^. In Fig. |9}3, we show the f^-behaviour of rj'' and r^^ 
given T. We deduce at the stability regions: 



r'^ccY = 



= 1.026(4)j,(4),„,(6)„,,.(l),^ , 



(34) 



where the errors coming from other parameters than ss are neg- 
ligible. This implies: 

Mai^ - Mh;^ = 94(27) MeV , (35) 

where we have used Mh^^ ^ 3.58 GeV from Eq. ( |3TI ) and the 
experimental value of M^^ . We show in Fig. [TO]the analogous 
analysis for the bottom channel. We deduce: 



r'\bby = 



Mn- 



= 1.0050(3),-,(10)„,,(4),(10)„, , (36) 



where the error is again mainly due to (ss), the others being 
negligible. This implies: 



Mn- -Mz- =50(15) MeV, 



(37) 



where we have used Mhj^ ^ 9.96 GeV using our prediction 
in the previous section. This resuh agrees with the potential 
model one of about 60 MeV given in ifioll . Again like in the 
case of spin 1/2 baryons, the SU(3) mass-differences appears 
to behave like the inverse of the heavy quark masses, which 
can be inspected from the QCD expressions of the two-point 
correlator One can also observe that the mass-splittings are 
almost the same for the spin 1/2 and spin 3/2 baryons. 



8. The £lftc/Sic mass ratio 



expressions of these corrections are: 



ImF 



\ I pert 



128 71^ fi 



6 £1 [(my - linlm])] 



+6 Z2mlf- - Abc^lt^ + (5m^ - 4m^) t 



-mt + Stnlml + Imt 



m,6(l -Hfe^) 



2^ 



32 nt^ 



Abc[P- + (ml + ml)t 



1 .007 
1 .0068 
1.0066 
1 . 0064 
1 .0062 



Figure 11: ili,c/Si,c : h behaviour of r^'^ibc): dot-dasiied line (red); r^ibc): 



continuous line (green), and r'^'^(bcy. dashed line (blue) for 
T = 0.8 GeV-2, nic = 1.26 GeV and nth = 4.22 GeV. 



1.00775 
1 . 0075 
1.00725 
1 .007 
1.00675 
1 . 0065 
1.00625 
1 .00 



50 GeV^ 




Figure 12: £ljc/Hjc a): r-behaviour of r[''(/5c): dot-dashed line (red) and 
rff(bc): continuous line (green) for k = -0.05, tc 



GeV. b): ^.-behaviour of rf{bc) and rf(bc) for r = 0.9 GeV"- and k 



50 GeV^ and = 1.26 
-0.05. 



The E(bc) and the Q.(bc) spin 1/2 baryons can be described by 
the corresponding currents: 



Je,,, = ea/iA [(clCysd/i:) + k(clCdp)y5] b,i. 



(38) 



where d, s are light quark fields, c, b are heavy quark fields and 
^ is a priori an arbitrary mixing parameter The expression of 
the corresponding two-point correlator has been obtained in the 
chiral limit m^ - ms = Oby Refs. We have checked 

these expressions which we complete here by adding the OTj- 
corrections for the PT and quark condensate contributions. The 



ImF: 



2 I perl 



ImF! 



2 \ss 



where: 



-2(ml-mlf] + —[(ml+my 

^hc 



-2(ml + ml)(ml - m^ft + (ml - m])"] 



(39) 



3msmcmt,(l + k^) 



64 n^t 

2mlml] - 2 X2 {ml - mfj t 
-Ahc [s + ml + mlj 
msmi,R(l - k^) 



2J:,[(ml + ml)t- 



Abe [t - ml + ml) + 



-^[mlf + (m^ + mlm] - 2mfj t 

+ {^1 - '«?)^] 



(40) 



Amlm^ 



(t - ml - m^y- 



1I/2 /. 2 2\ 



■£1 = :^ los 



1 +v 
1 - V 



-£2 = los 



(ml + ml)t + (ml - m^XA^^ - ml + m]) 



2mbmct 



(41) 



Like in previous sections, we study the different ratios of mo- 
ments in Figs. [TT]and[T2l As one can see in Fig. [TTk , rf(bc) 
and r^(bc) are quite stable in k and present common solutions 
for: 



k = +0.05 



(42) 



inside the range given in Eq. (|5]), while r^^^(bc) does not inter- 
sect with the other DRSR. The t and tc behaviours given in Fig. 
[T2h,b are also very stable from which we deduce the DRSR: 

r'''(bc) = ^ = 1.006(0.2),,(1.4)„Xl)me , (43) 

where the errors coming from other parameters are negligible. 
This implies: 



Mn,,-MH,„, =41(7) MeV, 



(44) 
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where we have used the QSSR central value M^,^^ ^ 6.86 GeV 
in Eq. Q- The size of the mass-splitting can be compared 
with the potential model prediction about (70-89) MeV given 
in lEIil. 

Table 2: QSSR predictions for the doubly heavy baryons mass ratios and sphttings, which 
we compare with the Potential Model (PM) range of results in 1 10, 15]. The PM prediction 
for the spin 3/2 is an average with the one for spin 1/2. The mass inputs are in GeV and 
the mass-splittings are in MeV. 



Mass ratios Mass inputs Mass plittings PM 





1.0167(19) 


--rr 


= 3.52ri3] 






= 59(7) 


70-93 


^y^bb = 


1.0019(3) 




= 9.94r9] 




- ^bb 


= 19(3) 


30-38 




1.0260(70) 


^cc 


= 3.52ri3] 


Cl'' 


— life 


= 92(24) 


90-102 


^bb/^bb = 


1.0049(13) 


^hb 


= 9.94r9] 


sib], 


- Sfci, 


= 49(13) 


60-73 


SI' /H* = 


1.0260(75) 


— ff 


= 3.58 •> 




~cc 


= 94(27) 


91-100 


'•'■bb'^bb 


1.0050(15) 


"66 


= 9.96 *' 


%b 


^ ^hb 


= 50(15) 


60-72 


Clbc/^bc = 


1.0060(17) 




= 6.86riO] 


Slbc 


- ^bc 


= 41(7) 


70-89 



*' We have combined your results for the mass-splittings with the experi- 
mental value of Mz^.^ and with the central value of M^^^^ in Eq. 0. 



9. Conclusions 

Our different results are summarized in Table |2] and agree 
in most cases with the potential model predictions given in 
MM- 

- The mass-splittings between the spin 3/2 and 1/2 baryons, de- 
rived in Eqs. (l22l i and ( l23l l is essentially due to the radiative 
corrections in our approach and seems to behave like 1 /Mq. 

- For the SU(3) mass-splittings, our results, derived in Eqs. dSTl i 
and (l33l l for the spin 1/2 and in Eqs. (l35T l and (l37T i for the spin 
3/2, indicate that the splittings due to the SU(3) breaking are 
almost independent on the spin of the heavy baryons but ap- 
proximately behave like 1/Mq. These mass-behaviours can be 
qualitatively understood from the QCD expressions of the cor- 
responding correlators where the leading mass corrections be- 
have like inJinQ. 

- Finally, we obtain, in Eq. ( l44l l. the SU(3) mass-splittings be- 
tween the Q.(bcs) and E(bcd), which is about 1/2 of the potential 
model prediction. 

Our previous predictions can improved by including radiative 
corrections to the SU(3) breaking terms and can be tested, in 
the near future, at Tevatron and LHCb. 
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